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! Abstract 

I In this work we identify the manifold of sohtary waves arising in a three- 

■ component scalar field model using the Bogomol'nyi arrangement of the energy 

^— I , functional. A rich structure of topological and non-topological kinks exists in 

^ I the different sub-models contained in the theory. 

o ■ 
^ ■ 

5 : 1 Introduction 

^ ■ The search for solitary waves is an ongoing topic in both Mathematics and Physics 

+^ ' because this kind of quasi-soliton plays an important role in a huge number of branches 

. of non-linear science. In Field Theory, they usually appear in models that support spon- 

. . ! taneous symmetry breaking, the most prominent examples being kinks/domain walls, 

I vortices, and monopoles Starting with theories that involve a high number of fields, 

^ I the usual procedure followed to investigate the existence of solitary waves -topological 

' defects- is to obtain an effective scalar field theory, imposing severe restrictions on the 

original theory. In most cases, one is compelled to pursue an effective theory that will 
correspond to a single scalar field model, where the existence of topological defects can 
be checked easily. The reason for this is the good-understanding of this kind of system; 
as paradigmatic examples, the well-known kink and soliton in the one- dimensional 0^ 
and sine-Gordon models should be noted. Both kinds of solitary wall can be thought 
of as thick walls, the topological defects in a three-dimensional perspective. Neverthe- 
less, the general framework is that the effective theory depends on several scalar fields 
and thus the truncation may involve an important loss of information concerning the 
presence of topological defects and the structure of spontaneous symmetry breaking. 
It is therefore desirable to investigate the general properties of domain wall solutions 
in a multi-scalar field theory. 
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In (l+l)-dimensional field theory, solitary waves are non-singular solutions of the 
non-linear coupled field equations of finite energy such that their energy density has 
a spacetime dependence of the form: e{x,t) = e{x — vt), where v is the velocity 
of propagation. In relativistic theories, Lorentz (or Galilean) invariance provides all 
the kink solutions from the purely static ones. The search for finite-energy static 
solutions in one- dimensional field theories is tantamount to the search for finite action 
trajectories in a natural dynamical system where, the x-coordinate plays the role of 
time; the field components transmute to positions in the configuration space, and 
the field theoretical density energy becomes minus the mechanical potential energy. 
No wonder the difficulties involved in finding kinks in multi-component scalar field 
theories: one faces multi-dimensional mechanical systems where integrability is not 
ensured. 

At a very early stage in the (pre-)history of the subject, a (l+l)-dimensional field 
theoretical model with two real scalar fields became relevant. Montonen and Sarker- 
Trullinger-Bishop proposed the deformation of the 0(2)-linear sigma model with a 
potential energy density of ?7(0i, ^2) = |[(0i + 02 ~ + '^^^i]; 0- clear 
that the zeroes of the potential are two points and hence the hunt for kinks started 
immediately^. Using a trial orbit method in the associated two-dimensional mechanical 
system, Rajaraman identified two different topological kinks joining the two vacua 
of the system that live on a straight line and half-ellipses respectively. Only one 
component of the scalar field is non-zero in the first case, but the two-components 
differ from zero in the second kind of solution; for this reason, these solitary waves are 
referred to as TKl (straight line) and TK2/TK2* (upper/lower half-ellipse) kinks in the 
literature that appeared later. Rajaraman also found one kink associated to a closed 
trajectory starting from and ending at the same point of the the vacuum orbit. Magyari 
and Thomas ^ realized that the mechanical system associated with the MSTB model 
is integrable -there is a second invariant in involution with the mechanical energy- 
and used this fact to show that there exists a whole family of two-component non- 
topological kinks (NTK2), all of them degenerated in energy with Rajaraman's NTK2 
kink; explicit kink form factors were only described by numerical methods. 

The main breakthrough in analytically finding all the solitary waves of the MSTB 
model emerged in Reference Ito discovered that the mechanical problem was not 
only integrable but that it was Hamilton- Jacobi separable by using elliptic coordinates. 
In this setting, he showed the analytic formulas for the kink orbits and the kink form 
factors, unveiling the mathematical reasons for the previously observed striking kink 
sum rule. Immediately, the stability of this degenerate kink family was questioned; 
application of the Morse index theorem solved this problem in jS] . A parallel with the 
Morse theory of geodesies was established somewhat later in Reference [0]. Thus, a 
clear connection arose between solitary waves, their stability, and dynamical systems. 
In Reference [7], several of us showed that the MSTB model is not unique in this 
respect; two (l-l-l)-dimensional field theoretical models with two real scalar fields - 
referred to as model A and B in that paper- have manifolds of solitary waves with 

^We shall refer to the zeroes of the potential as vacua throughout the paper, anticipating their 
role in the quantization of this classical field theory. Also, because these two points are related by the 
internal symmetry group Z2 (81 Z2 generated by (f)i and (f)2 -~4>2, we shall sometimes refer to 

this set as the vacuum orbit. 



2 



similar structures. To find the analytic expression for the kinks of model A, we were 
prompted to solve an integrable dynamical system classified as Liouville Type I, see 
[H]. The system belongs to the same class as that found in the MSTB model -the 
two-dimensional Garnier system [H]- but there are three differences: (a) the potential 
energy density is a polynomial of sixth order in the fields (instead of fourth); (b) the 
vacuum orbit has five points (instead of 2), and (c) there are many more stable kinks 
than in the MSTB model. Model B is characterized by a fourth-order potential energy 
density in the two scalar fields. The main feature, however, is the need to solve a 
Liouville integrable system of Type III, i.e. Hamilton- Jacobi separable in parabolic 
coordinates. The vacuum orbit has four points and there are manifolds of stable and 
unstable kinks. 

In recent years, all this work has proved to be fruitful in the framework of super- 
symmetric theories. In the dimensional reduction of a generalized Wess-Zumino model 
with two chiral superfields, Bazeia-Nascimento-Ribeiro- Toledo (henceforth referred to 
as the BNRT model) [10] found one one-component topological kink (TKl) and one 
two-component topological kink (TK2). In this case, the vacuum orbit has four points 
and the potential energy density is a polynomial in the fields of order four. Under- 
standing the BNRT model as a deformation of model B, some of us discovered the 
whole manifold of kink orbits ^l] . There is kink degeneracy, also found slightly earlier 
by Shifman and Voloshin in one of the topological sectors and, for two critical 
values of the coupling constant, analytic formulas for the kink form factors are avail- 
able. One of them corresponds exactly to model B; the other one leads to a Liouville 
system of Type IV, Hamilton- Jacobi separable in Cartesian coordinates. Interesting 
consequences have been translated to the dynamics of intersecting branes [TB]. How 
thick walls grow from one-component kinks is well known. Composite kinks give rise 
to a non-trivial low energy dynamics for intersecting walls as geodesic motion in the 
kink moduli space (the space of the integration constants with a metric inherited from 
the field theoretic kinetic energy). Another supersymmetric model that shows a rich 
pattern of kink solutions is the Wess-Zumino model itself. The BPS kink states of this 
M = 2 supersymmetric (1 -|- 1)D model with a complex scalar field and holomorphic 
superpotential were discovered by Vafa et al. in [T3j. In (TH], two of us studied this 
system from the point of view of the real-analytic structure. The vacuum orbit having 
been identified, the flow between the vacuum points was determined as the gradient 
of the real(imaginary) part of the superpotential. Thus, kink orbits are identified with 
real algebraic curves. 

Here, we continue to struggle with the extension of these studies to field theoretical 
models with three real scalar fields. In some of us explored the generalization 
of the MSTB model. The solution of the three-dimensional Garnier system using 
three-dimensional Jacobi coordinates revealed the existence of an extremely complex 
variety of kinks. Nevertheless, the structure of the kink manifold and its stability was 
completely unraveled in J7] . The main goal of the present paper is to identify the kink 
manifold arising in a family of three-component relativistic field models with a vacuum 
manifold that contains several elements or points. This family can be interpreted as the 
natural generalization of the generalized MSTB model studied in ^Hl QZ| in the sense 
of Stackel-type systems. The most interesting feature of this generalization is that the 
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number of elements in the vacuum manifold depends on the range of relative values 
of the coupling constants. Therefore, we can find different submodels of our system, 
which have a very rich structure of kink manifolds. When the energy density of these 
kinks or solitary waves is studied we find that several families of these solutions are 
degenerate, which allows us to claim that some kink families indeed consist of more 
basic kinks, such that their energy density displays several lumps associated with the 
basic kinks. In our model, we are able to find solutions with two, three or four lumps. 

The organization of the paper is as follows: In Section 2 we introduce the model, 
writing the expressions in Stackel form and describing the different spontaneous symme- 
try-breaking scenarios. Section 3 is divided into four sub-sections. In 3.1, we identify 
first-order differential equations satisfied by the kink solutions, reproducing the Bo- 
gomol'nyi procedure in this context. Sub-section 3.2 contains the resolution of these 
equations. In 3.3, we determine the regions where the solutions live and, finally, in 
Sub-section 3.4, some general comments about the determination of the stability of 
the kink solutions are offered. In Section 4 we describe the behaviour of solitary wave 
families in one of the regimes of the model, at the same time discussing their stabil- 
ity properties. Finally, in Section 5 we address some points concerning the different 
extensions of the model. 



2 The model 



We focus our attention on the search for kink solutions arising in three-component 
scalar field models in a (H-l) Minkowskian space-time, whose dynamics is governed 
by the action functional 



S[<P] 



1 ^ 
2^ 



where we use Einstein's convention for Greek indices with the usual metric rjn = —r]22 = 
1, V12 = V21 = 0, and where is a smooth non-negative function that depends on the 
three-component scalar field = 025 03)- We use natural units, hence c = 1, and 
we shall henceforth denote x^ = t and x^ = x. The Euler-Lagrange equations in this 
case are written as the following system of second-order partial differential equations 

920,; ^20,; dU _ 



dx"^ 



^1, V2, 



1,2,3. 



Kinks are finite-energy solutions of ((H), such that the time dependence is dictated by 
the Lorentz invariance: (j)K{t,x) = <P{^f^)-, and they can be interpreted as extremals 
of the positive semi-definite energy functional 



£\ 



dxe(x) 



dx 



IE 



dx dx 



+ f/(01,02,03) 



(2) 



which maintains this functional finite: S[(f)] < +00, see PJ. Therefore, solitary waves 
must comply with the asymptotic conditions 

d(f) 



lim I 

X— >itoo 



e M 



h) 



lim 

x^itoo dx 







(3) 
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where Ai is the set of zeroes or absolute minima of the potential term -that is, Ai = 
{(01,02,03) £ K^/^(0i, 02, 03) = 0}- which are usually referred to as vacua of the 
theory because the elements of A4 play this role in the corresponding quantum theory. 

The usual procedure for tackling the search for kinks in this kind of theory is to 
interpret Q as the action functional of a mechanical system in which we think of the 
variable x as "time"; as the coordinates of a unit-mass point particle, and V = —U 
as the potential function. From this point of view, are merely equations of motion 
in the new system. In reference j^, the authors deal with the model involving the 
potential function 



and show that the mechanical analogue is not only completely integrable but also 
Hamilton- Jacobi separable by using a system of three-dimensional elliptic coordinates. 
In [T7], the stability properties of kinks are analyzed and a new approach to search 
for kinks based on the Bogomol'nyi decomposition are given in the above system. 
The authors prove the equivalence between the Hamilton- Jacobi equation and the 
Bogomol'nyi approach. The potential function (jH) has two zeroes, = (—1, 0, 0), v'^ = 
(1, 0, 0). Therefore, the kinks in this model can be classified into topological and non- 
topological kinks according to whether the solution connects two different vacua (open 
orbits) or the solution departs and arrives at a vacuum (closed orbits). 

The search for new integrable models is not an easy task. In this sense, we would 
remark the following quotation from Jacobi in his "Vorlesungen liber Dynamik" , which 
allows us to see the issue from a different perspective: "The main difficulty in integrat- 
ing given differential equations is to introduce suitable variables which cannot be found 
by a general rule. Therefore, we must go in the opposite direction and, after find- 
ing some remarkable substitution, look for problems to which it could be successfully 
applied". 

The goal of this paper is to generalize the above model, focusing our attention on 
models with a greater-than-two number of elements in A4, such that we can find a 
more sophisticated symmetry-breaking scenario and a richer plethora of solitary waves 
than before. 

Using the same notation as in the reference ^Hl, we now introduce a system of 
Jacobi elliptic coordinates A = (Ai,A2,A3), with constants a| = 1 — a|, (t| = 1 — 
and 1, which is defined as: 



[/(0i, 02, 03) = l{<Pl + <Pl + <Pl- ir + ]-al<Pl + lal<Pl 



(4) 





-1 




'I 



-1 
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3 



cr: 



in which the range of the coordinates is: 



- oo < Ai < 0-3 < A2 < 0"2 < A3 < 1 . 



(6) 
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It should be noted that this coordinate transformation is invariant under the group 
G = Zf generated by 0a ^ (-l)^"''0a; b = 1,2, 3. 

Invoking (jSJ, the energy functional can be written as 



dXj dXj 
dx dx 



+ f/(Ai,A2,A3)^ , (7) 



where the metric coefficients QniX) = — ttt — ,w/^^-I/, — have been introduced. 

3 

Here, we set /j(A) = JJ^ (Aj — A^). 

In the new variables, the potential (j^ is written as: 

and their zeroes v~ y f + are mapped to one point v = {X\, Ag, A3) = (0, (t|, a"^) in the 
elliptic space because of the above-mentioned invariance. 

In order to generalize expression (jHl), we introduce the following potential function 

U{X; a^) = J2 U.iX; a') = ' " , (9) 

i=i i=i •^^^ ^ 

which becomes a polynomial function of eighth degree in the original fields. Notice 
that we have added a new factor (Aj — a^)^ to each of the summands in (jH)). Thus, Q 
introduces new degenerate vacua in Ai, which for fixed a| and depend upon the 
value of the coupling constant = 1 — a^. Therefore, new scenarios of spontaneous 
symmetry-breaking and a richer kink manifold arise in this model. Taking into account 
the range © for the elliptic coordinates and formula 0, we can observe that the new 
structure of the set A4 depends on the relative values between the constant and the 
fixed constants ct|, cxg and 1. For instance, for > 1 the new factor (Aj — a^)^ does 
not vanish for any value of Aj and therefore Ai has the same structure as that in model 
dHI). However, for G (cr|, ct|) we find new vacua located at the points A = (ct|, a^, ct|) 
and A = (0, a^, a-|). 

We shall now introduce different scenarios for our model depending on the value of 
the constant a^. We shall distinguish the number of vacua in each case. 

• Regime El: As mentioned above, for G (1, 00) there exists only one vacuum 
in the elliptic space, minimizing the potential function: A''^ = (0, 03^, (f2^)- We 
have a similar situation if the constant takes the discrete values 0, a| or 0^2. 
For this reason we define the set Lq = {0,a|,a2} U (l,oo), taking into account 
that if G Lq our model only has a vacuum, A^^, in the elliptic space. In the 
Cartesian space, the vacuum manifold Ai can be regarded as the orbit generated 
by the action of the group G/Hi over the vacuum vi, where ifi = 1 x Z2 x Z2 
is the group that leaves the coordinates of vi invariant. There are therefore two 
vacua in the Cartesian space A4o = {(j)" =(±1,0,0)}. 
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The kink solutions in this model display the same behaviour as those of the model 
studied in jTHj , although the explicit expression of the equations of motion is more 
complicated because we have a polynomial of degree eight in the original fields. 
Owing to this similarity, we shall not deal with this regime in our study. 

• Regime E2: We now consider the range G Li = {OjO'D for the coupling 
constant. In this regime, new zeroes of the potential arise on the plane Ai = a^, 
in the elliptic space that corresponds to the ellipsoid + -/^-2 + ./^-a = 1 

in the Cartesian space. In fact, two vacua arise in the elliptic space, A*"^ = 
(0, (Ts^, (72^) and A^ = (a^, 03^, (T2^), both invariant under the subgroup H2 = Hi. 
Correspondingly, there are four vacua in the Cartesian space that correspond to 
the orbit ULi('^/^i)^i- Therefore, we have Mi = {^""^ = (±1,0, 0),^^^ = 
(±a,0,0)}, as depicted in Figure 1. 

It is interesting to remark that the range of values G (—00, 0) is formally 
analogous to that in which G Li, interchanging the roles of the ellipsoids 
Ai = and Ai = in the previous reasoning. We shall therefore focus our 
attention on the range of values Li. 




Figure 1: Vacuum manifold in the Cartesian and elliptic spaces: Regime E2. Fi, F2 and F3 
stand for the foci of the ellipsoid; B and C are the extremes of the minor semi-axis, and A 
represents the umbilical points. 



• Regime HI: In this case, G L2 = {a\,a'^. From the values of and the 
range of the elliptic coordinates, the zeroes of the potential term ^ arise on 
the plane A2 = a^, which is equivalent to the hyperboloid of one sheet + 
-i-n = 1 + in the Cartesian space. We find three vacua located at the 

points A''^ = (0,a|,a|), A''^ = (a-|,Q;^,a-|) and A"^ = {Q,a^,al). The vacua 
Vi and V2 remain invariant under ifi, whereas fs is invariant under i^a = 1 x 
Z2 X 1. There are eight vacua in the Cartesian space corresponding to the orbit 
[£=ii.G I Hi)vi, with coordinates M2 = {(j)"' = (±1, 0, 0), ^''^ = (±a, 0, 0), ^''^ = 
(±^, 0, ±^-\/a^ — (t|)}, as shown in Figure 2. In section 4, for the sake of 
clarity we shall focus on this regime in order to describe in detail a particular 
kink manifold of the model instead of discussing it in each single regime. 

• Regime H2: This case is characterized by G L3 = (ag, 1). Applying the 
same reasoning as before, we find that new vacua arise on the plane A3 = a^; 
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Figure 2: Vacuum manifold in the Cartesian and elliptic spaces: Regime HI. 



that is, the hyperboloid of two sheets 



1 + 



+ 



in Cartesian 



coordinates. In particular, the potential has four minima: A^^ = (0,(T3,(J2;, 



\Vi _ (r> ^2 ^2^ 

= (^^2^ ^2^ ^2^. ^ (0,a^,a2)^ and A''^ = (0,a-|,a2). The vacuum is 
invariant under H^^ = 1 x 1 x Z2, and the Cartesian vacuum manifold is the orbit 
\At=iiG/Hi)vi; namely, 

M, = {0^^ = (±1,O,O),0'^^ 



(±«,o,o). 



0-3 



0-3 



^2 



cr 



0"2 0"2 



ai,0 




B, 
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Figure 3: Vacuum manifold in the Cartesian and elliptic spaces: Regime H2. 



Regime H2': In this case, the coupling constant is set equal to unity. In 
the internal elliptic space we can read the minima as A^'^ = (0,(T3,(t|), A*'^ = 



(erg, (t|, 1), A^^ = (0, ct|, 1) and A"-* = (0, (T3, 1), which correspond to eight minima 



h^3 



in the Cartesian space: M3 = lim M3 = {0"' = (±1, 0, 0), ^''^ = (0,0,0), 
(O,O,±cT3),0^* = (O,±a2,O)}. 

In this latter case, the plane A3 = 1 is introduced into the elliptic space. Unlike 
the previously introduced planes, this is no longer a regular one, and this can 
be readily seen in the degeneracy exhibited by the H2 vacuum manifold at the 
limit — 1; this singular plane corresponds to the plane 0i = in Cartesian 
coordinates. Regarding the kink manifold, this is basically the same as that of 
the H2 model, except that the kink solutions existing on the two sheets of the 
hyperboloid and in between them now degenerate into kink solutions on the plane 
01 = 0. This situation is the 3D analogue of model A in [Zj. 
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3 First-order equations and Kink Manifolds 



3.1 The superpotential and the Bogomol'nyi arrangement 

We notice that the potential determines a Stackel system jH]. Therefore, the 
Hamilton- Jacobi equation of the mechanical analogue is separable using the system 
of Jacobi elliptic coordinates. However, here we shall make use of the Bogomol'nyi 
arrangement in order to obtain the kink manifold of our model. The two procedures 
are equivalent (see ^'^^ the second one allows us to identify the supersymmetric 
extension of our field theory, given that if the energy functional ((Tj) can be written as 

' ' '''^(^' (10) 



dx 



for some function W{Xi), then the underlying field theory has a supersymmetric exten- 
sion in which the function W plays the role of superpotential in the supersymmetric 
field theory, see [IH]. Therefore, the superpotential W must comply with 



dXj ' ''''' ' \dX2J ' \dXs 

Plugging the expression of the potential function ^ and the metric coefficients into 
the above equation, we have 

^ Af (A. - ar UU(^^ - -4(A. - 1) UU(^^ - (dW- ' 

h m h m Ua, 

which can be solved easily by the ansatz W = Wi{Xi) + W2{X2) + Wsi^X^). The three 
resulting decoupled ordinary differential equations 



dW,y _X^iK-a% ^--12 3 



'_ X^{K-a ^^^ 
dXi J 4(1 - A 

lead us to the expression of the superpotential function W 

3 3 

wi^M)^x) = J2wt{x,) = -J2i-'^y'P2{x^)V^^^ , A = o,i, 

i=l 1=1 

where P2(Ai) = 2d + dXi - 3A^, with d = {5a^ - 4). 

Extremal trajectories for the energy functional ()10|) arise if the following system of 
first-order differential equations 



^_^^,.^Hk-a^)(\-ol)(\-<rt)^^^ (11) 

Ji\X) 



9 



where /3j = 0, 1 and i = 1, 2, 3 is satisfied, because the squared terms in ^TU^ are always 
positive and the last one is a constant. Due to the indeterminacy of the signs f32 
and Ps, (|TT|l constitutes eight systems of ordinary differential equations. Nevertheless, 
this set of systems is easier to solve than second-order (Euler-Lagrange) equations. In 
order to obtain a complete kink solution we have to join solutions from the first-order 
differential equations with different choices of the signs (—1)^' in different intervals 
covering the real line. The reason for this is that the first-order differential equations 
inherit the information of the second-order equations defined piecewise. Assuming 
that we search for continuous and differentiable solutions, the sequence of signs (—1)^' 
corresponding to the different pieces that constitutes a solution is prescribed. In section 
4 we shall illustrate this approach in several cases. From ()10|) it is readily seen that 
the energy of a solitary wave, solution of with only one piece, depends only on the 
topological charge of the solution. In this case, it is said that the Bogomol'nyi bound 
is saturated. However, if the orbit A is given by A = U^^^A-^, where J is the number of 
pieces of A and A-^ stands for the j^^ piece, we have: 

pieces of A j «=i 

= i:(^^''^^(^final)-^^''^^(^Litial))' (12) 
j 

where represents the values of the Pi parameters for the j^^ piece of the solution. 



3.2 Solutions via quadratures 

In order to solve system (fTT|l . we rewrite it in the form: 

d\i dx 

(-i)/^»2v/r^nl=i(A. - 9) ^ m ' 

ined c 

of these equations gives 



2 = 1,2,3, (13) 



where we have defined c = (a^, (t|, cr|, af) and cr| = for future convenience. The sum 



y 4 =0 (14) 



Multiplying each side of (fT^ by Aj and summing over i, we obtain: 

y =0. (15) 

tr(-i)^»2v/r^n;=i(A.-c,) 

Also, multiplying (fT!?|l by Xf and summing again over i we reach the equation that 
establishes the dependence of the kink components on x 

y = dx. (16) 



1=1 
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We shall now determine the kink orbits and the form factor by invoking ()14|) . (fT3j) . 
and (|T6|l . Integration of the first two equations, 



i=l 

3 

E 



-1) 



ft 



2 



XjdXi 



VI- Ai ]\j=i{\i-Cj) 



72 



73 



leads us to the expression of the generic kink orbits: 



,272 



.273 



n 



n 



(-1)^1 



n 



^T^+V^ 



n 



(-1)^3 



(17) 



n 











\/i-A2+iyi-cj 





3 

n 



{-l)'53c 



1— c. 



where Fj{c) = a/1 — Cj Y[i=i ly^ji^j ~ '^O; 72 and 73 are arbitrary real constants 
that specify a particular kink orbit. 

The integration of ()16|) 



(-l)^-> 



VI- Aini=i(Ai-cj) 



7i + a; 



gives us the form factor of the kink: 



,2(71+2:) 



n 

i=i 





^ ' J 3 




(-l)/52c2 






Vl-Ai-^l-Cj 


' ' n 






Vl-A3--\/l-Cj 


F,{c) 











(19) 

7i being an integration constant associated with the translational invariance of the 
system. Expressions f|T7j) . (|T8|) and (|T9|) provide us with the whole manifold of sohtary 
waves. 



3.3 Frontiers and barriers. Basic kinks 

We shall now prove that the set of solitary waves is confined to living in a bounded 
region of the internal space, which in fact corresponds to a parallelepiped in the elliptic 
space. For the sake of clarity, we shall restrict our study to the range G L, where 
L = [J^^i Li is the set in which the kink manifold is richest, see Section 2. This include 
the regimes E2, HI, and H2. Squaring the first equation in (fT^ . and defining the 
generalized momentum tti = gii(\)^, we have: 



^,2N2 



-vr? - M(^^i^ = 0. (20) 
2 ^ 8 1-Ai 
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Equation ()20|) can be regarded as that governing the motion of a particle moving 
under the influence of the potential function 



2\2 



W(Ai 



(1-Ai) 



3 



, — OO < Ai < (j^ 
oo , 0"? < Ai < oo 



The function has at least one minimum in Ai = and a second one in Ai = 
if G Li. Furthermore, the function W(Ai) goes to — oo as Ai tends to — oo. 
Thus the bounded motion can only occur in the interval [0,ct|]. This, combined with 
the boundary conditions, leads us to the conclusion that the kink solutions lie in the 
parallelepiped P^iO) = [0,0-3] ^ [^3' ^2] ^ [^i'l]- 

There is still more information that can be extracted following this procedure, owing 
to the appearance of a second minimum. Let us first fix a value in L, and let us 
set G Li for some i that depends on a^. Squaring the i^^ equation of the system 
and defining the generalized momentum tTj = gii(X)^, we arrive at a similar 
one-dimensional dynamics: 

1 , A?(A, - 5^)^ 
2"' - 8(1 -A.) = °- 

Accordingly the corresponding potential function is now defined by U{\i) if z = 1 
and 



W(A, 



A? (A. - a'? 
8(1 -A,) 



min{Li} < Aj < max{Lj} 



00 , Ai ^ Li 



for i = 2,3. The minimum Aj = now separates the bounded motion of the one- 
dimensional system into two intervals - the Aj G L^ = [min{Lj},a^] interval and the 
Aj G Lf = [a^, max{Lj}] interval -, and into the trivial motion Xi = L^ = . This, 
together with the asymptotic conditions, leads us to conclude that, besides living in 
P3(0), the kink solutions lie entirely in the sets 

P3(0)-'0'+ = {A G P3(0) with A, gL-'0'+}. 

This decomposition of the parallelepiped -P3(0) is, for the case we shall study in 
detail, regime ifl, as follows (see Figure 2): 

P-m = hiO)- U P3(0)° U P3(0)+ = [0, al] X [^3^, a'] x [al 1] U 
U [0, al] X {a'} X [al 1] U [0, ^3^] x [a^ aj] x [aj, 1] . 

The parallelepipeds P3(0)~ and P3(0)^ contain families of solutions that depend 
on two and three parameters, whereas the plane -^3(0)*^ only contains two-parametric 
solutions. 

Thus, introduction of the factor (A, — a^)^ into the potential function U{X) leads us 
(within our range of study) to a new confinement of kink solutions in the parallelepiped 
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PaiO). The generic kink solutions divide into two sectors and, in addition to this, a 
new kind of two-parametric solutions arises: those satisfying Aj = a^. Consequently, 
the kink manifold can be decomposed as follows: 

C = CruC°UC+, (22) 

where Cj~ represent the class of kink solutions with Aj < a^, Aj = and Aj > 
respectively. 

3.4 Stability 

In this sub-section we discuss how to determine the stability properties of the kink 
solutions. For the whole variety of kink solutions in this system, it is not possible to 
solve Ai, A2 and A3 in terms of elementary functions of x. Therefore, it is not possible 
to explicitly write out the Hessian operator for any kink in the model and, hence, the 
stability properties cannot be studied through analysis of its spectrum. 

To determine the stability of the solutions, we use instead the arguments developed 
in Ref. based on the Jacobi fields along kink solutions. Although the treatment 
depicted in that paper is for a deformed Sigma 0(3) model, the extension to this 
model can be readily carried out. Following this procedure, a rule establishing the 
stability (instability) of the solutions is obtained: each solution crossing either the 
edge F1F3 = {ct|, cr|, A3} or the edge AF2 = {Ai, o"!, ct|} becomes an unstable solution, 
since these two edges constitute lines of conjugate points of each vacuum of the theory. 

The key point is that the superpotential function is not different iable over either 
of these two edges and, consequently, the energy of the kink (fT^ is not a topological 
quantity since it depends on the value of the superpotential at the crossing points. 

In what follows, and bearing this remark in mind, we shall only mention the char- 
acter of each of the kinks described. 

4 Description of the Kink manifold in the HI regime 

The description of the kink manifold in the different regimes arising in our model 
is a long and tedious task. We shall therefore focus our attention on a particular 
example: the HI regime. Nevertheless, this case will suffice to illustrate the general 
features that also arise in other regimes of our model. We shall now describe the 
behaviour of the kinks that arise in the HI regime of our model. We can find basic 
kinks, similar to the solutions TKl and TK2 in MSTB model, that are placed on the 
edges of the characteristic parallelepiped in the elliptic space (see figures 4,5 and 6). 
These solutions are the simplest kinks in our model and they consist of a single lump, 
such that they can be interpreted as an extended particle. We shall show that the 
kink manifold includes other kink solutions involving several lumps associated with 
the basic kinks. 

We recall some remarkable points of the HI regime from the previous sections: The 
number of minima is three in the "elliptic" space, and eight in the Cartesian one (see 

Figure 2): M2 = {^^ = (±1,0,0), = {±a,0,0),<j)^^ = (±^, 0, ±f - ^1)}- 
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The ellipsoid -E = 0? + ^ + % = l (that is, A = (0, A2,A3)), the one-sheet hyper- 

boloid H = + — ^23^ = 1 (o^ -^2 = a^), and the planes ^2,3 = are 

distinguished surfaces in the internal space. In 3.3, we have proved that all the topo- 
logical solutions are confined within the above-mentioned ellipsoid E. From this point 
of view, these surfaces play the role of separatrices among three-parameter families of 
solutions, as proved above. These solutions are associated with finite values of the in- 
tegration constants, 7^. It is usual in the literature [12] to refer to this class of solutions 
as generic solutions. On the other hand, these surfaces also contain the trajectories of 
two-parameter families of solitary waves, which correspond to asymptotic values of the 
constants 7^. Accordingly, they are called non-generic solutions. 

Finally, we describe the kink manifold in these cases. We can distinguish: A, Non- 
generic, two-parametric families, and B, Generic, three-parametric families of solitary 
waves: 



A Two-Parametric families of solutions: 



Al Solutions on the ellipsoid E. 

The potential term Ui vanishes on this surface. Accordingly, the superpo- 
tential function is: 



3 3 
$^W,'''(A.) = -5^(-lf^P2(A.)v/r^ , A = 0,1. 



W^(''^'*^(A2,A3) = 

1=2 1=2 

The orbit of these solutions is given by 

(-l)^2ej 

^2,2 A Vl^2- ^^-^ 



3 

n 



Co- 



C4 



3 

n 



VI - A3 + ^/l - Cj 



72 being an arbitrary real constant. We have two kind of solutions: 

i) T^'^^ : Stable topological solutions that connect the minima and 

after having crossed the plane 0i = 0. 
ii) N""^ : Unstable non-topological solutions that join the minimum with 
itself. The trajectory of these solutions starts from ^3, reaches the plane 
01 = 0, and -after crossing the umbilical point A- returns to the same 
point f^; see Figure 4. 



The energy of these solutions can easily be calculated by integrating dW 
along their respective orbits: 



/ dW= f + / 



dW^ + 2 / dW. 



a 



— a 



(To 



2,f-. 
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Figure 4: Solitary waves on E in the Cartesian (left) and elliptic (middle) spaces. Energy 
density of a kink of the family T^''"^ (right). 



dW^ + 4 



dW^ 



15 



(-2(T2P2(^2)-2a^2(a')) 



(To 



— + a" 

5 / V 5 



A2 Solutions on the plane 03 = 0. 

In this case, the terms U2 and Ui of the potential vanish, but not simulta- 
neously. The former vanishes over A2 = and the latter over Ai = (t|. 
Because of this, two superpotential functions appear, and hence two sys- 
tems of differential equations must be involved in order to determine this 
solution. Nevertheless, we can synthesize W as follows: 

W(''-^^nAfc,A3) = -^$^(-l)^'P2(A.)v/r^ , A = 0,1, 



=fc,3 



where = 1 for A2 = o"!, and k = 2 for Ai = (t|. The equations of the orbit 
on the plane A^ = (Xg are: 



,272 



n 

i7^3 



^1^- 



(_l)/3fc(c._c3) 



n 



(-l)^3(c.-C3) 



Again we have two kinds of solutions: 

i) T^^i''"^ : Unstable topological solutions linking the vacua and v"^. 
These solutions leave f^, intersect the axis 02 and the segment 
consecutively, and finally arrive at w^, as depicted in Figure 5. 
A^l : Unstable non-topological solutions connecting v"^. The solutions 

0, cross the focus F2, and return to 



go from V , intersect the axis 02 
the initial point v"^. 



15 



^3 S e(x) 




X 



Figure 5: Solitary waves on (/>3 = in the Cartesian (left) and elliptic (middle) spaces. 
Energy density of a kink of these families (right). 



The computation of the energies is as follows: 



(73 J 



— a — 



— a 



A3 Solutions on the plane 



0, see Figure 6. 



Now, the terms and U2 of the potential vanish over A3 = (t| and A2 = o"2; 
respectively. The two superpotential functions that appear can be synthe- 
sized in a similar way: 



0,1, 



i=l,k 



where A; = 2 for A3 = (j| and k = 3 for A2 = The equations of the orbit 
on the plane A^ = are: 



,272 



n 

37^2 



(_l)/3l(c.-c2) 



n 

jV2 



(-l)/3fc(e^.-e3) 



We now have three classes of solutions: 

i) T^l'^'^ : Stable topological solutions that join the minima and f^, as 

can be observed in Figure 6. 
ii) TJ|: Unstable topological solutions that connect the point with the 
minimum, which is its reflection by the transformation 03 —(ps, pre- 
viously crossing the focus F3. 
iii) TJ2'^3- Unstable topological solutions that link the points v"^ and v^. 
In this case, the solutions depart from f^, and finally arrive at f ^ after 
intersecting the axis 03. 
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(Jj ^3 e(x) 




X 



Figure 6: Solitary waves on </>2 = in the Cartesian (left) and elliptic (middle) spaces. 
Energy density of a concrete kink of the family (right). 



The energies for these solutions are: 

'1 



0-2 J 



m 
e[t:i 



5"" 



a 



a 



providing a simple kink energy sum rule: 2E\T^ 
remaining energy is: 



a 



a 



V\,V2\ 



(72 J 



The 



f2,^'3l 
0-2 -I 



(To 



-ai] -2 



a 



— a 



— a 



In figure 6(right), we have depicted the energy density e{x) of a member of 
the family . We notice that the kinks of this family consist of three basic 
lumps. 

A4 Solutions on the hyperboloid. 

The term U2 vanishes over A2 = and hence the superpotential function 
is: 

W-(''^'*HAi,A3) = ^$^(-l)''^P2(A.)v/r^ , A = 0,1. 
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The equation of the orbit is: 



4 = 1,3 



,272 



4 

n 

i=2 



(_l)/3l(c.-ci) 



4 

n 

j=2 



1 - Cj 



(_l)/33(c,-ci) 



In this case, only one family is found. 



* Tj^' ^ : The trajectories of these stable solutions connect the points v 
and f ^, previously intersecting the plane 0i = 0, as is shown in Figure 7. 
Notice that the energy density in this case comprises two basic lumps. 

The energy is: 



«2,f31 



(To 



(T. 



a 
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B Three-Parametric families of solutions. We find three kinds of solutions: 

Bl Solutions located inside the ellipsoid and outside the hyperboloid, see Figure 
8: 

i) T^^'^^ : Stable topological solutions that join f ^ and f^. The solutions 

emerge from v^, later cross the plane 0i = 0, and finally arrive at v"^. 
ii) T'"'-* : Unstable topological solutions, which start from a minimum v^, 
consecutively cross the planes (pi = and 03 = 0, intersecting the FiF^ 
edge, and finally arrive at v^. Notice that the energy density in this 
case comprises four basic lumps. 




Figure 8: Generic solitary waves in the Cartesian (left) and elliptic (middle) spaces. Energy 
density of a kink of the family T""^ (right). 



a 
5 



Their energies are: 
B2 Solutions located inside the hyperboloid: 



— a — 



a 



— a 



a 

5 
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i) T^2,V3 . xhese are unstable solutions. They leave v"^, cross the plane 
01 = 0, later intersect the hyperbola AF2, cross the plane 0i = again, 
and finally arrive at the point v^; see Figure 9. 




Figure 9: Generic solitary waves in the Cartesian (left) and elliptic (right) spaces. 
The energy in this case is: 




To complete the previous energy calculations, the kink energy sum rules satisfied 
by the generic solutions are offered: 



See Sub-section 3.1 of Reference ^Hj for an explanation of the origin of these 
rules in a simpler setting. We stress that the decomposition of the kink energy 
density in several lumps is due to the kink energy sum rules. 

Finally, as an example we depict the kink form factor (Fig. 10 and Fig. 11) for the 
two unstable generic solutions. 

5 Further Comments 

It is possible to generalize this kind of model in two senses; we enlarge the internal 
space with N scalar fields and we include a greater number of coupling constants af. 
1. To study the generalization of this kind of system to dimensions, it is first 
necessary to introduce A^-dimensional Jacobi elliptic coordinates. An appropriate ex- 
planation of these can be seen in ^E] . The potential function we propose for the system 
is as follows: 



i=l i=l 
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Figure 11: Form factors in the Cartesian space for the T'"^ and T'"^''"^ solutions 
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where the couphng constants together with the coordinates satisfies the chain 



— oo < Ai < (T^ < A2 < . . . < Ajv_i < (J2 < Aat < 1 = cr^. 

The denominator is /i(A) = nj^il-^i ~ -^j); ^^"^ ^ ^ positive constant. The 
function f/(A; a^) is positive semi-definite and presents a number of zeroes, depending 
on . The most interesting kink manifold appears when G Lj , i = 1, . . . , iV, and 
becomes richer as i increases, the Lj intervals being the trivial generalization of those 
appearing in the three-dimensional potential. 

We shall now briefiy study the vacuum manifold in all the N different cases at once. 
Let us set o? such that G Lj for some j between 1 and A^. To find a zero of the 
function f/(A; o?\ we must make every term f/i(A; o?^ vanish. To visualize the process, 
we shall seek help from the following graphic 

9 k ^ ■■■ I ■■■ I I 



j + 1 N-j 



Each circle in the A^ block represents a value that A^ can take to make the term 
f/fc(A; a^) null. Each value appearing in the vacuum coordinates will be represented by 
a full circle, and hence each vacuum in the elliptic space is represented by N full circles. 
To fill the N — j circles to the right of a^, there is only one possibility, as seen in the 
figure, but to fill the remaining j circles we have a number of different ways equal to 
the number of permutations of j -|- 1 elements, j of them being repeated. Therefore, 
we have Pj^i~^^ = J + 1 zeroes of the U{X; a^) function, j of them being on the plane 
Xj = o? . To figure out the number of corresponding Cartesian vacua, we only need to 
take into account the multiplicity of each elliptic vacuum. By doing this, we conclude 
that by introducing a regular plane A-,- = cP' there are V = 4 -|- (j — 1)2^^^ Cartesian 
vacua. The kink manifold thus decomposes into disconnected sectors 0- 

2. The second generalization considers not only one parameter, o? ^ but several of 
them. The generalized potential is constructed as follows. 

Let us consider numbers nj = 0, 1, 2, . . . with i = 1, . . . , A^, and let us define [rix -|- 
. . . + n7v) different parameters a\p such that for each 7^ 0, a\j G Li and j = 1, . . . , nj. 
We can therefore construct the A^-dimensional potential: 

A'" 1 -'^ \'2,(\ -2\/\ -2\ 

t/ ^ ^ [/.(A. 4) = 1 ^ - ■ - n - 4)^ ■ (23) 

i=l i=l •^^^ ' j=i 

The case in which 'Yl!i=\"^^ ~ corresponds to the deformed 0{N) linear sigma 
model [in] and the case J2iLi = 1; with = 3, is precisely the model studied in the 
previous sections. 

As YliLi increases, the vacuum manifold becomes more and more abundant owing 
to the appearance of an increasing number of roots in the potential. An easy way to 
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account for the vacuum manifold V is through the corresponding generahzation of the 
previous graphic 



'^i I I '^i I Vi I I '^j I '^j+i I I '^N-i I ^: 



)oo...«( 




In this picture {rii + ... + 11^) additional circles appear, a-holes for short, since 
every afj is easily seen to be a root of the Ui term in ()23p. Computation of the number 
of vacua now proves to be an easy task given that, as before, each vacuum point is 
represented by N filled circles. It happens that the number of vacua -including vi, 
which corresponds to all the a-holes emptied- is given by: 



Card(V) = 1 + ^ AT, 

q=l 



N 

M 



where Ng is the number of vacuum points with q filled a-holes, which can be calculated 
readily using combinatorial techniques. 

Regarding the kink manifold, and looking at the corresponding first-order equations, 
for each afj we can deduce a confinement of the solutions in PsiO) similar to that 
obtained in section 3. Therefore, a number of 2("i+---+"-^) subsets of -P3(0) that host 
general kink solutions appear. 

The purpose of this construction is now clear. Recalling the stability criterion and 
the confinement of the solutions due to the factors (Aj — af^)^, we can isolate the edges 
-^1-^3 = {^3; (^3, ^3} and AF2 = {Ai, (t|, Proceeding in this way, we can find subsets 
of the configuration space in which only stable solutions emerge. 

Acknowledgements 

The authors warmly acknowledge many discussions with J. Mateos Guilarte. A.A.I, 
thanks the Secretaria de Estado de Educacion y Universidades of Spain for financial 
support. 

References 

[1] R. Rajaraman, Solitons and instantons. An introduction to solitons and instantons in 
quantum Held theory, North-Holland Publishing Co., Amsterdam, 1982 

[2] C. Montonen, Nucl. Phys. B112 (1976) 349-357 

S. Sarker, S.E. Trullinger and A.R. Bishop, Phys. Lett. A59 (1976) 255-258 
K.R. Subbaswamy and S.E. Trullinger, Physica D2 (1981) 379-388 

[3] E. Magyari and H. Thomas, Phys. Lett. AlOO (1984) 11-14 

[4] H. Ito, Phys. Lett. A112 (1985) 119-123 



22 



[5] H. Ito and H. Tasaki, Phys. Lett. A113 (1985) 179-182 

[6] J. Matcos Guilarte, Lett. Math. Phys. 14 (1987) 169-176 
J. Mateos Guilarte, Ann. Phys. 188 (1988) 307-346 

[7] A. Alonso Izquierdo, M.A. Gonzalez Leon and J. Mateos Guilarte, J. Phys. A31 (1998) 
209-229 

[8] A. Perelomov, "Integrable Systems of Classical Mechanics and Lie Algebras", Birkhauser, 
Boston MA., 1990 

[9] R. Garnier, Ren. Circ. Mat. Palermo 43 (1919) 155-191 

[10] D. Bazeia, J.R.S. Nascimento, R.F. Ribeiro and D. Toledo, J. Phys. A30 (1997) 8157- 
8166 

D. Bazeia, H. Boschi-Filho and F.A. Brito, JHEP 9904 (1999) 028 

[11] A. Alonso Izquierdo, M.A. Gonzalez Leon and J. Mateos Guilarte, Phys. Rev. D65 
(2002) 085012 

[12] M. Shifman and M. Voloshin, Phys. Rev. D 57 (1998) 2590 

[13] A. Alonso Izquierdo, M.A. Gonzalez Leon, J. Mateos Guilarte and M. de la Torre 
Mayado, Phys. Rev. D66 (2002) 105022 

[14] S. Cecotti, C. Vafa, Comm. Math. Phys. 158 (1993) 569-644 

R Fendley, S.D. Mathur, C. Vafa, and N.P. Warner, Phys. Lett. B243 (1990) 257-264 

[15] A. Alonso Izquierdo, M.A. Gonzalez Leon and J. Mateos Guilarte, Phys. Lett. B480 
(2000) 373-380 

[16] A. Alonso Izquierdo, M.A. Gonzalez Leon and J. Mateos Guilarte, Nonlinearity, 13 
(2000) 1137-1169 

[17] A. Alonso Izquierdo, M.A. Gonzalez Leon and J. Mateos Guilarte, Nonlinearity, 15 
(2002) 1097-1125 

[18] A. Alonso Izquierdo, M.A. Gonzalez Leon and J. Mateos Guilarte and M. de la Torre 
Mayado, Ann. Phys. 308 (2003) 664-691 



23 



